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SUMMARY 


Six months of data (Sept. 1968 through Feb. 1969) from the wide-f ield-of- 
view low resolution infrared radiometers (LRIR) on the Environmental Science 
Services Administration (ESSA) 7 satellite were analyzed. Earth-emitted and 
Earth-reflected irradiances were computed at satellite altitude using data from 
a new in-flight calibration technique. Flux densities and albedos were computed 
for the top of the Earth’s atmosphere (30 km). Monthly averages of these quan- 
tities over 10°-latitude zones, each hemisphere, and the globe are presented for 
each month analyzed, and global distributions are presented for typical months. 
Emitted flux densities are generally lower and albedos higher than those of pre- 
vious studies. This may be due, in part, to the fact that the ESSA 7 satellite 
was in a 3-p.m. Sun-synchronous orbit and some of the comparison data were 
obtained from satellites in 12-noon Sun- synchronous orbits. ESSA 7 detectors 
seem to smooth out spatial flux density variations more than scanning radiome- 
ters or wide-f ield-of-view fixed-plate detectors. Significant longitudinal and 
latitudinal variations of emitted flux density and albedo were identified in the 
tropics in a zone extending about ±25° in latitude. 


INTRODUCTION 

Solar energy plays the dominant role in generating the Earth’s weather and 
climate. Our capability to study this role took a giant step forward with the 
advent of artificial Earth satellites, since they provide excellent platforms 
from which to monitor the amounts of solar energy which the Earth reflects, 
absorbs, and reemits. Starting with Explorer 7 in 1960, instruments for measur- 
ing components of the Earth’s radiation budget have been flown on about half the 
United States meteorological satellites. Earth radiation budget results from 
the first generation of these meteorological satellites (1960 through 1965) were 
reported in reference 1 . 

Included in the second generation of meteorological satellites, which began 
in 1966, was the ESSA series. Radiation measurements made with the low resolu- 
tion infrared radiometers (LRIR) on the ESSA satellites are an important subset 
of the radiation measurements of the late 1960's. Descriptions of techniques 
for analyzing these data and results from very limited analyses of the radiation 
data from ESSA 3 are given in references 2 and 3, but no previous paper has 
directly reported results from radiation analyses of the ESSA 7 satellite data. 
One possible reason that the ESSA data have been largely ignored was a detector 
degradation problem which was caused by continued exposure to solar radiation. 

In references 2 and 3, techniques are described for overcoming this problem 
using in-flight data, and an improved technique for flight calibrating the detec- 
tors was reported in reference 4. The new technique requires neither the assump- 
tions concerning the behavior of emitted flux density at spacecraft ingression 
and egression that are required by the techniques of references 2 and 3 nor the 
assumptions relative to net flux density required in reference 3. 



This paper describes results of an analysis of ESSA 7 radiation data using 
the in-flight calibration technique of reference 4. Irradiances are computed at 
the satellite, and flux densities and albedos are determined for the top of the 
Earth's atmosphere. Flux densities and albedos are presented for each month of 
data analyzed, and global distributions of these parameters are shown for typi- 
cal months. ESSA 7 results are discussed relative to results from several pre- 
vious radiation budget analyses, and some results are compared directly to other 
results. 
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SYMBOLS 

albedo, ratio of reflected to incident solar flux density 

area of ith 1 0°-latitude zone, m 2 

detector housing constant 

ir radiance, W/m 2 

geometric shape factor 

solar flux density at mean distance from Earth to Sun, W/m 2 

measurement (satellite) altitude, km 

reference altitude (top of atmosphere) , 30 km 

thermal mass coefficient of detector, J/m 2, K 

distance from Earth to Sun, normalized by mean distance 

radiant flux density or exitance, W/m 2 

radius of Earth, 6371 km 

temperature, K 

absorptance of detector 

detector response factor 

spherical coordinate defined in figure 18 

angle between spacecraft spin axis and Sun, deg 

emittance of detector 

spherical coordinate defined in figure 1 8 
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H max unrestricted view angle of Earth horizon from detector 

0 latitude 

X longitude 

K angle incoming ray makes with normal to detector surface 

o Stefan-Boltzmann constant, 5.66961 x 1 0“ 8 W/m 2 *K 4 

$ Sun zenith angle at subsatellite point, deg 

do) solid angle at detector subtended by increment of area 

Subscripts: 

A ascent 

av average 

b black-coated 

D descent 

d detector 

e Earth emitted, longwave 

G globe or entire Earth surface 

h detector housing 

n net 

r Earth reflected, shortwave 

s direct solar 

w white-coated 

Other notations: 

cos* <|> function (defined in fig. 5) for computing albedo 
Bar over quantity denotes monthly average. 

Dot over quantity denotes derivation with respect to time. 
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SATELLITE ORBIT AND DATA SET CHARACTERISTICS 


Orbit 

The ESSA 7 meteorological satellite was launched by NASA on August 16, 1968. 
The orbit was Sun synchronous and nearly circular with an average altitude of 
1450 km, a period of 114.7 min, and an inclination of 102°. The satellite 
crossed the equator during ascent (going northward) at a mean local time of 
3 p.m. The average local time is shown in figure 1 as a function of satellite 
latitude. Measurements were made every 30 sec, or about every 1.6° of Earth 
central angle. At an altitude of 1450 km, the spacecraft is in direct sunlight 
at subsatellite-point Sun zenith angles up to about 123°, whereas the terminator 
at the Earth's surface occurs at a Sun zenith angle of 90°. Ingression and 
egression as used herein will refer to the spacecraft's entry into and exit from 
direct sunlight. This always occurs on the descent portion of the orbit so that 
during ascent the spacecraft is continuously exposed to direct solar radiation. 


Data Collection 

Measurements were made almost continuously from August 16, 1968, to 
July 19, 1969. The data set obtained from the Environmental Data Service 
of the National Oceanic and Atmospheric Administration (NOAA) covered the 
period from Sept. 1, 1968, to May 31, 1969. Data for the second two weeks 
of December were missing from this set, and sampling for the first two weeks 
of December was not as complete as for the other months. Due to anomalies 
in the data of March, April, and May, 1969, results reported herein are for 
the period September 1968 through February 1969 only. 

Because the orbit is Sun synchronous, the local, or Sun time, at which the 
spacecraft crosses a given latitude, on ascent or descent, is essentially con- 
stant. This means that the average emitted flux density at a given latitude is 
based on measurements made at only two local times (one on ascent, one on 
descent) . Average reflected flux densities and albedo at the same latitudes 
are based on ascent measurements, only, and thus include measurements at only 
one local time. In references 1 and 5, variations in emitted flux densities 
with local time were found. A significant effect on albedo due to local time 
was also found in reference 1 . Many of the other results discussed in this 
paper are based on data sets obtained with Sun- synchronous orbits. The data 
of reference 6 were obtained with a 3-p.m. orbit, those of reference 7 were 
obtained in a 12-noon orbit, and some of the data of references 1 and 8 were 
obtained with satellites having various Sun- synchronous orbits. 


Data Editing 

Four different data editing schemes were used to eliminate bad or uncertain 
data. The first scheme, described in reference 4, eliminated points when the 
measurements were obviously inconsistent with the orbit ephemeris or when there 
were inconsistencies in the ephemeris itself. For instance, measurements might 
indicate that the spacecraft was in direct sunlight while the ephemeris showed 
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the spacecraft was in darkness, or the spacecraft position would be inconsis- 
tent with the Sun zenith angle. 

The second editing scheme dealt with data in the transition regions where 
rapid heating and cooling of the detectors at ingression and egression produced 
uncertainties in the measurements. Data in these regions were rejected when 
the Sun zenith angles were in the range from 123° to 127°. The zones affected 
were always on the descent portion of the orbit, and for the affected zones, 
the distributions may not be centered in the zone. 

In the third editing scheme, data were eliminated when the computed emit- 
ted ir radiance was found to lie outside a given envelope. For each latitude 
zone, an envelope was given which was based on a computed value of the esti- 
mated mean and variations about the mean. This scheme takes into account the 
significant pole- to-equa tor gradient in the mean emitted flux density as well 
as changes in the variation about the mean from zone to zone. 

The fourth editing scheme eliminated data when the magnitude of the com- 
puted irradiance changed more than a specified amount in a specified time. 

This scheme can reject spurious data spikes where the magnitude of the irradi- 
ance could be within the envelope specified in the third editing scheme. 


Spatial Sampling 

Spatial sampling, both in number of measurements and in distribution, is 
believed to be adequate for each of the six-month periods analyzed. Contours 
showing the spatial distribution of the edited sample for October 1968 (which 
is typical) are presented in figure 2 for both the ascent and descent portions 
of the orbit. These contours are based on the number of samples in regions of 
10° latitude by 10° longitude. The distribution is seen to be quite uniform 
during ascent. The less uniform distribution during descent is due mainly to 
the elimination of data at spacecraft ingression and egression. 


DESCRIPTION OF DETECTORS AND MEASUREMENTS 

ESSA 7 was a spin-stabilized satellite with its spin axis normal to the 
orbit plane. Each detector consisted of a pair of flat-plate, bolometer-type 
radiometers wired in electrical series and mounted opposite each other on the 
outer circumference of the spacecraft. The temperatures of the detectors T^ 
and of the housings which support the detectors T^ are the basic measured 
signals. A more complete discussion of the ESSA 7 detectors can be found in 
references 2 and 4. Each radiometer alternately views the Earth and space, 
and the integrated signal from the two radiometers is treated as an average 
measurement over the rotation period of about 6.5 sec. The instrument system 
was made up of two detectors. Each detector has different surface property 
(called white or black) to discriminate longwave- and shortwave-radiation com- 
ponents. Each detector is in radiation equilibrium. That is, 

Ed = M d (1 ) 
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where E d , the irradiance absorbed by the detector from external sources, is 
given by: 


^ = (a/s ) e E e + (a/e) r E r + (a/e) s E s (2) 


and M d , the exitance at the detector, is calculated, using the measured temper- 
atures of the detector and its housing, by the equation 


4 ( 4 4\ 

M d = a T d + ac(T d - T h j + K*T d 


(3) 


The equation for the irradiance absorbed by the detector (eq. (2)) can be sim- 
plified in two ways. 

First, a detector's response to Earth-reflected radiation and to direct 
solar radiation will be the same since these two radiation sources have essen- 
tially the same spectral characteristics. That is. 


(a/e) s = (a/e ) r 


Second, we will assume that we can independently compute a value for the 
direct solar irradiance E s . Now if we let M d ^ w be the calculated value of 
radiant exitance for the white-coated detector (eq. (3)) and if we make use of 
the equilibrium relationship (eq. (1)), we have 


M d,w “ (° t / e )e,w®e + ( a / G )r,w®r + ( a / e )r,w^s (4a) 

Similarly for the black-coated detector, we have 


^d,b “ (&/ c )e,b®e + ( a / G )r,b-^r + ( a / e )r,b®s (4b) 

Equations (4a) and (4b) can be solved simultaneously for the Earth-emitted 
and Earth-reflected components of irradiance, E e and E r . 


Md,b ( a / G ) r ,w M d,w( ot / £ ) r ,b 

g ~ _ . _ 

( a / G ) r f w ( a / £ ) e f b ” ( a / e ) r ,b ( a / £ ) e,w 


(5a) 
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e,b M d,b( a / e ) e,w 

g — , , 

r (a/e) r ,w( a / £ ) e r b ~ (a/e ) r , b ( a / e ) 


(5b) 


E s 

e r w 


The section which follows describes how the detector characteristics 
required to compute emitted and reflected irradiances at the satellite were 
determined by equations (5a) and (5b). 


DETECTOR CHARACTERISTICS 

Coefficients Required for Computations of Detector Exitance 

Values for the detector housing constants C required to calculate the 
detector exitance by equation (3) were determined in preflight laboratory cali- 
bration tests. The calibration procedure, which consisted of testing a flight 
version of the detector, is described in reference 4. The resultant detector 
housing constants were found to be 


C w 0.075 

C b 0.054 


The detector thermal mass coefficients K' were calculated using flight 
data obtained during spacecraft ingression and egression. The calculations are 
described in reference 4, and the resulting thermal mass coefficients are plot- 
ted in figure 3 as a function of detector temperature. Exitance due to changes 
in detector temperature (thermal mass term in eq. (3)) was small outside the 
regions of ingression and egression because temperature changes were gradual. 
Abrupt temperature changes inside these regions, however, resulted in sharp 
changes in exitance, and some measurements had to be rejected. 


Detector Response to Longwave Radiation 

Determination of absorptance and emittance of the two (black and white) 
detectors for emitted radiation was made for several temperatures and included 
spectral absorptance data for the two detector coatings. The spectral absorp- 
tance curves and the values of absorptance and emittance as functions of tem- 
perature are given in reference 4. One value of longwave absorptance was 
selected for each detector. Each value corresponded to the temperature when 
the Earth was assumed to be an equivalent blackbody emitter. Two values of 
longwave emittance were selected for each detector. These corresponded to the 
detector's response when in and out of direct solar radiation. The table on 
the following page gives the values of detector absorptance a e and the two 
corresponding values of the detector response factor (ot/e) e for each of the 
two detectors. The process for determining the ingression and egression points 
in the orbit, which indicate when to make the changeover in detector response 
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factors, is described in reference 4. The selection was based on the Sun zenith 
angle at the subsatellite point. The temperature change of the detectors was 
used as a backup indicator. 


Detector 

ot „ 
e 

(a/e ) e 

Direct sunlight 

No direct sunlight 

White-coated 

Black-coated 

0.80 

.94 

1 .053 
1 .000 

1 .081 
1 .021 


When the spacecraft is not in direct sunlight, the last two terms of equa- 
tions (4a) and (4b) are zero because only Earth-emitted radiation is being mea- 
sured. For this condition, the ratio of the exitances in equations (4a) and (4b) 
should give the ratios of the detector's longwave response factors. When ratios 
of the computed exitances were plotted for the no-direct-sunlight condition, the 
values were found to be nearly constant and equal to the ratios as determined 
from the selected values of (0t/£) e . This gives confidence in the technique for 
selecting values of longwave absorptance and emittance. 


Detector Response to Shortwave Radiation 

Based on flight test (refs. 2 and 9) , the type of coatings used on the 
ESSA 7 detectors has a history of degrading. Ground tests at NASA Langley 
Research Center have shown the degradation was primarily a result of exposure 
to ultraviolet radiation. This degradation produces changes in the response 
of the detectors to shortwave radiation, and the white coatings were found to 
be much more vulnerable to degradation than the black coatings. 

In the current analysis, the Sun is used as a calibration source to deter- 
mine values of the shortwave-detector response factors. Just before the space- 
craft exits from direct sunlight (egression) , a detector measures longwave 
and direct solar radiation. After entry into total darkness, the detector 
measures longwave radiation only. If it is assumed that longwave radiation 
is constant for these two conditions, the difference in computed exitances 
for the white detector (using eq. (4a)) is 


AM d/W = <aA) r ,w E s 


and similarly for the black detector 


^d,b ~ ( a / £ )r,b E s 
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The procedure can also be followed at ingression, and thus, values of 
(a/e) r ^ w and (<*/£) r,b can be inferred for the detectors twice during each 
orbit. A modified version of the preceding method was developed (ref. 4) which 
does not require the assumption that the longwave radiation remains constant 
for the period of the calculation. The values of (a/e) r determined by this 
method for the two detectors are shown in figure 4. Values for the white- 
coated detectors are seen to vary about 25 percent over the data period, whereas 
values for the black-coated detector vary much less. These changes, although 
quite significant, are much smaller than those reported in reference 6 for 
the ITOS 1 instrument which changed its response by about 40 percent in a 
two-month period. 


COMPUTATIONS 

Irradiances 

Values of emitted and reflected irradiances at the satellite were computed 
by equations (5a) and (5b) for each measurement in the edited data set. Values 
of solar irradiance E s required in equation (5b) were computed by 


E 


s 



sin y 
TT 


The value of the mean solar flux density H s is 1 353 W/m 2 . Values of H s 
used in previous studies have ranged from 1 353 W/m 2 (ref. 6) to 1 395 W/m 2 
(ref. 1). A fairly recent experiment, reference 10, estimated a value for H s 
of 1367 W/m 2 . 

The geometric shape factor sin y/ir relates the component of solar irradi- 
ance measured by the detector to the incoming solar flux density. A derivation 
of this shape factor is given in the appendix. The value of y varied between 
42° and 48° during the data period, and there was a corresponding change in 
the latitude — local-hour relationship. These variations reflect the departure 
from a truly Sun- synchronous orbit and the seasonal variation in the Sun’s 
declination. 


Flux Densities 

Irradiances at the satellite were reduced to flux densities at an altitude 
of 30 km (top of the atmosphere) above the Earth (assumed radius = 6371 km) by 
the relationships 


M e = E e /F 

(6a) 

M r = E r /F 

(6b) 
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where F is the geometric shape factor. Values for the geometric shape factor 
were calculated from the function derived in the appendix. It was assumed in 
the derivation that the Earth's radiance field was Lambertian. It was shown 
in references 11 and 12 that this assumption is valid for estimating emitted 
flux densities. Reflected radiation is generally not diffuse at Sun zenith 
angles greater than about 30°. However, most previous estimates of reflected 
flux densities have been made using the diffuse assumption. The calculated 
value of the shape factor was 0.2433 at a mean satellite altitude of 1450 km 
and varied about ±2 percent with orbit altitude. These values were found to be 
in very close agreement with values obtained using the methods of references 13 
and 1 4 . 


Average Values of Flux Density and Albedo 

Monthly averages of emitted flux density were derived from values obtained 
from equation (6a) for 10°-latitude zones and for 10°-latitude-by-10°-longitude 
regions for both the ascent and descent portions of the orbit. The average 
local times corresponding to these averages of emitted flux densities can be 
obtained from figure 1 . For example, the average local times for ascent and 
descent in the 0° to 10° zone are approximately 3 p.m. and 3 a.m. Monthly 
averages of zonal and regional emitted flux densities based on the combined 
data set were obtained by 


^ = 


M e, A + M e,D 
2 


(7) 


This averaging method precluded biasing the monthly averages for the combined 
data set due to unequal sampling between ascent and descent. 

Monthly averages of zonal and regional albedos were computed by 


(M r )av 

A = 

Hg (1 /R) 2 cos* $ 


( 8 ) 


where ( M r^av the fl ux density in the zone or region obtained by averaging 
the reflected flux densities of equation (6b) . The denominator is the average 
incident solar flux density corresponding to the measurements. It represents 
the average reflected flux density that would be measured for an albedo of 
unity. This incident solar flux density differs from that which would have 
been obtained using cos <J) because the Sun zenith varies over the detector 
field of view. The difference is greatest at Sun zenith angles near 90°, as 
can be seen in figure 5. This is because a large portion of the detector 
field of view is not illuminated by the Sun. Average albedos were not 
reported when the average zenith angle was greater than 90°. 
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Monthly averages of zonal and regional reflected flux densities were 
computed by 

M r = M S A (9) 

where A is obtained from equation (8) . The monthly average of the incident 
solar flux density M s for the zone or region is obtained by integrating 
incident solar flux density over the entire zone or region. 

Monthly averages of zonal and regional net flux density were computed by 


M n = M s - (M e + M r ) = M s (l - A) - M e 

Monthly averages of hemispherical and global emitted, reflected, and net flux 
densities were obtained by area weighting the monthly averages on zonal scales. 
For example, the monthly average of global emitted flux density is 


Me,G 


I 
1 =1 




I 
1 =1 


a i 


Monthly averages of hemispherical or global albedo are the ratios of the monthly 
averages of reflected flux to the monthly averages of incident solar flux. For 
instance, a monthly average of global albedo is 


18 


Ag " 


^ ( M r,i) a i 
i=l 

18 

^ <M s ,i)ai 


i=l 


n 


Three-month averages of hemispherical and global flux densities are obtained 
by taking the average of the three separate monthly averages. The three-month 
average of hemispherical or global albedo is the ratio of the three-month 
average of reflected flux to the three-month average of incident solar flux. 

This aver-aging technique prevents biasing the three-month averages toward 
months with larger data sets. 


RESULTS AND DISCUSSION 
Global and Hemispherical Averages 

ESSA 7 .- Monthly and seasonal averages of hemispherical and global flux 
densities and albedo are given in table I for the ESSA 7 data period. Global 
emitted flux densities have a maximum and a minimum in October and January f 
respectively, with corresponding variations from the six-month average of about 
±2 percent. Month- to-month changes in global emitted flux density during the 
winter months are due almost entirely to changes in emitted flux density in the 
Northern Hemisphere. Global albedo has a maximum, which occurs in November. 

The maximum variation from the average global albedo was less than 4 percent 
during the data period. This variation is due to changes in the albedo of the 
Southern Hemisphere. Averages of hemispherical net flux density are dominated 
by hemispherical solar flux density. For instance, albedo in the Northern 
Hemisphere has only a slight variation during the fall months, and emitted flux 
densities vary only about 4 W/m 2 for the three months. Net flux density, how- 
ever, changes from 21 W/m 2 in September to -57 W/m 2 in November due to the sea- 
sonal decrease in hemispherical solar flux density. 

Comparisons with other results.- Presented in table II are seasonal aver- 
ages of global flux densities and albedos from ESSA 7 and three other sources 
(refs. 1, 8, and 15). Some global results from references 6 and 7 are also 
shown for comparison with appropriate ESSA 7 results from table I. Seasonal 
averages of the global emitted flux density from the reference sources are seen 
to be from 0 to 4 percent higher than those of ESSA 7, and there is up to 3 per- 
cent difference between averages within the three comparison sources. Combined 
averages for the two seasons are 225.5 W/m 2 for ESSA 7 and 231.8 W/m 2 for the 
reference sources. The corresponding blackbody temperatures are 251.1 K for 
ESSA 7 and 252.9 K for the references. The lower values obtained with ESSA 7 
data are not surprising because references 1 and 5 indicate that a Sun- 
synchronous orbit with a 3-p.m. equator crossing should result in relatively 
low values of emitted flux density. Emitted flux density data of Nimbus 3 
(ref. 7), which was in a 1 2- noon Sun- synchronous orbit, are higher than those 
of ESSA 7 by greater differences than those between ESSA 7 and the other three 
reference sources. The orbit of the NOAA 1 satellite (ref. 6) was nearly 
identical to that of ESSA 7, but the global average of emitted flux density 
from reference 6 is substantially greater than that of ESSA 7 for February. 

(See table I.) This may be due, in part, to assumptions necessary to determine 
degradation in the response of the NOAA 1 radiometers from in-flight data. The 
satellite was stabilized in three axes, and the detectors did not view the Sun 
directly. This meant that the Sun could not be used as an instrument calibra- 
tion source as was done with ESSA 7. 
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Averages of global albedo from ESSA 7 are seen to be higher in the fall 
period and lower in the winter period than those of references 1, 8, and 15. 
Averages from these reference sources also show an increase in albedo from 
fall to winter whereas ESSA 7 data shew a decrease. Very little difference 
is seen in albedo between the two periods shown from reference 7, and these 
values are both higher than the corresponding values for ESSA 7 in table I. 
Results of reference 6 show albedo values which are greater than those of 
all other sources discussed. Most of these differences are probably due to 
sampling differences imposed by the particular orbits, the majority of which 
were Sun synchronous. 

Fall averages of global net flux density are in good agreement with those 
of references 8 and 15, but the winter averages agree better with results from 
reference 1 . The large averages of emitted flux density and albedo obtained in 
reference 6 produced a loss in net flux density compared to a gain determined 
with ESSA 7 data. 


Averages in 1 0°-Latitude Zones 

Emitted flux density .- Monthly averages of zonal emitted flux density are 
given in figure 6. The largest month- to-month variations in emitted flux den- 
sity are seen for the fall months in the Southern Hemisphere, and the smallest 
variations are seen for the winter months in the Southern Hemisphere. Compar- 
isons of October and January averages for ESSA 7 with the results of reference 8 
are given in figure 7. ESSA 7 results indicate that latitudinal variations in 
the equatorial zones are not as sharp as those seen in the results of refer- 
ence 8, which included some data from scanning radiometers. The ESSA 7 detec- 
tors, which are rotating flat plates, generally smooth out the spatial varia- 
tions in flux density more than scanning radiometers or wide-f ield-of-view 
(WFOV) , fixed, horizontal flat plates. Averages of reference 8 are higher 
in the equatorial zones for both October and January, and agreement is reason- 
ably good at latitudes above about 35°. An exception is seen for the Southern 
Hemisphere in January where results of reference 8 are significantly higher than 
those of the current study. October and January averages of emitted flux den- 
sity are shown in figure 8 for the ascent and descent portions of the orbit. 
These generally show the ascent averages near the equator, which are taken near 
3 p.m. local time, to be higher than the descent averages, which are taken at 
3 a .m. 


Albedo .- Monthly averages of zonal albedo are shown in figure 9. Latitu- 
dinal variability in the equatorial region is less in fall than in winter, and 
the differences between fall months in these regions are very small. Between 
fall months, differences in albedo near the poles in both hemispheres follow a 
seasonal pattern. That is, as the Sun moves southward from September to Novem- 
ber, there is an increase in albedo in the Northern Hemisphere and a decrease 
in the Southern Hemisphere. In table I, note that the significant decrease in 
local albedo at high northern latitudes from December to February produces 
almost no change in hemispherical albedo. The opposite can be seen in the 
Southern Hemisphere. The accuracy of estimating albedo from WFOV measurements 
has always been a problem when the Sun is near the horizon, because the field 
of view of a WFOV detector includes some area which is illuminated by the Sun 
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and some which is in darkness. A parameter estimation technique described in 
reference 16 offers some promise of improvement in estimating reflected flux 
densities and albedos at high zenith angles. As with emitted flux density, dif- 
ferences in albedo between winter months near the South Pole are small. Compar- 
isons of October and January averages of albedo between ESSA 7 and results of 
reference 8 are given in figure 10. Agreement, particularly near the equator, 
is seen to be very good. 

Net flux density .- Monthly averages of zonal net flux density are presented 
in figure 11. The general shape of the monthly profiles and the month- to-month 
variations in net flux density are dominated by solar flux density. Maximum val- 
ues of net flux density lag maximum solar flux density by a few degrees of lati- 
tude, except for September where net and solar flux densities are coincident in 
latitude. The accuracy of the variation near the poles is questionable due to 
the problem mentioned earlier in estimating albedo near the poles. Also, the 
poles are not illuminated in some months. Comparisons of October and January 
averages of net flux density between ESSA 7 and reference 8 are shown in fig- 
ure 12. Agreement is very good in October and reasonably good in January. 


Global Distributions 

The monthly averages of flux densities and albedo in regions of 1 0° lati- 
tude by 10° longitude were used to generate global distributions (or maps) of 
radiation parameters. A cubic spline contour program was used with a spline 
tension factor of 2.5. The parametric cubic spline is generated by computing 
a polygonal contour through the defined grid and then generating the splined 
contour by using the polygonal contour length (normalized to one) as the inde- 
pendent variable. Contour slopes are matched when a contour closes. ESSA 7 
results are presented in figures 13 to 15. 

Contour maps of emitted flux density and albedo, presented in figures 16 
and 17, show comparisons of results from ESSA 7 (Sept. 1968) with results from 
Nimbus 6 (Aug. 1975). The Nimbus 6 values are previously unpublished results 
obtained by the authors. The data are from the wide-field-of-view detector 
of the Earth radiation budget (ERB) experiment on the Nimbus 6 satellite. (See 
ref. 17.) The contour technique is the same as that used on the ESSA 7 data. 

Emitted flux density .- Global maps of emitted flux density are shown in 
figure 13 for October and January. It can be seen that the ESSA 7 detectors 
identified significant structure in emitted flux density in a tropical zonal 
band extending about ±25° in latitude. Longitudinal gradients in this zonal 
band are, in some areas, as great as latitudinal gradients. For instance, there 
is a strong positive gradient extending westward from a low over western South 
America to a high over the Pacific Ocean. The influence of the structural band 
extends northward in October to include some of North America and the Mediterra- 
nean Sea, but the limit of the northward influence in January is below the tro- 
pic of Cancer. In this band, highs are generally centered over oceans and 
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lews over land f usually coinciding with minimum and maximum regions of cloudi- 
ness , respectively. Photographs taken by Vidicon cameras on the ESSA 7 satel- 
lite show a nearly cloud-free region in October corresponding to the indicated 
high in emitted flux density centered over southern Saudi Arabia, and the 
January photographs showed a buildup of cloudiness in that region. 

Comparisons of emitted flux density are given in figure 16. The struc- 
ture in the zonal band of influence is quite similar for the two data sets. 

The structure is more detailed for the Nimbus 6 data, however. Some of this 
may be due to the differences between September and August, but it is probably 
also due to the higher resolution of the flat-plate detector on Nimbus 6. The 
magnitudes of the Nimbus 6 data are also higher. This may be due to the differ- 
ence in local time when the two satellites cross the equator. The Nimbus 6 
orbit was Sun synchronous with a 12-noon northward equatorial crossing. 

Albedo .- Global maps of albedo for October and January are shorn in fig- 
ure 14. Comparison of figure 14 with figure 13 shows that areas of low albedo, 
in general, correspond to areas of high emitted flux density, and vice versa. 
This is more apparent in January when the longitudinal structure in the tropics 
is more clearly defined. It is interesting that the decrease in emitted flux 
density from October to January over southern Saudia Arabia does not have a 
corresponding increase in albedo, which is as sharp in detail. This may be due, 
in part, to the southerly shift of the band of longitudinal variability. Out- 
side this band, latitudinal variability overrides longitudinal variability. As 
with emitted flux density, the comparison maps of figure 17 show remarkable sim- 
ilarities between the September albedos from ESSA 7 and those for August from 
Nimbus 6. For instance, the highs in albedo over Africa are almost identical in 
shape for both data sets, and the contours on the east and west sides of South 
America are nearly the same for both data sets. 

Net flux density .- Global maps of net flux density for October and January 
are presented in figure 15. The dominant influence of solar flux density on 
the global distribution of net flux density can be seen in figure 15 where the 
width of the zone showing significant longitudinal variations is more narrow 
than for either emitted flux or albedo. 


CONCLUDING REMARKS 

The ESSA 7 data set is an important part of the Earth radiation measure- 
ments of the late 1960*s. For the six months of data analyzed, the time and 
spatial coverage was quite good. Results are enhanced by using the in-flight 
calibration technique developed in NASA CR-1 32623. The values of emitted flux 
density are lower and the values of albedo are higher than most corresponding 
values found in previous radiation budget studies. This may be due, in part, 
to the 3-p.m. Sun- synchronous orbit of ESSA 7. Latitudinal gradients dominate 
the spatial distributions of emitted flux density and albedo near the poles. 
Significant longitudinal gradients, however, were identified in a tropical zone 
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extending ±25° in latitude. Highs in emitted flux density occurred mostly over 
oceans and lows over land, while the opposite was found for albedo. Monthly 
changes in zonal and hemispherical net flux densities are dominated by changes 
in solar flux density. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 7, 1979 
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APPENDIX 


GEOMETRIC SHAPE FACTORS FOR THE ESSA 7 DETECTORS 
Upwelling Earth Radiation 

Let the radiance be Lambertian (diffuse) everywhere on a spherical surface 
above the Earth, and let the flux density M be a function of longitude X 
and latitude 0. An increment of irradiance at a satellite due to radiation 
from an increment of area on the surface can be written 


1 

dE = - M(X,0) cos £ dU) 
TT 


where ? is the angle which the incoming ray makes with the normal to the 
detector surface and dU) is the solid angle at the detector subtended by 
the increment of area. The total irradiance due to the Earth region within 
the detector field of view (FOV) is given by 


E = - C M(X,0) cos £ doo 

* J 

FOV 


If the flux density is assumed to be uniform over the detector field of view, 
the irradiance can be expressed as a product of the flux density and a geometric 
shape factor F. That is, 


E = MF 


where the geometric shape factor is given by 


■IS 


F = - \ cos £ dw 

FOV 


(AT) 
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APPENDIX 


Figure 18 illustrates the geometry of the ESSA 7 detector. The rotation angle 
3 of the detector was in a plane normal to the surface of the Earth. In this 
spherical coordinate system. 


and 


5 = n 


dw = cos ri dri d3 


Thus, equation (Al) can be written 


4 p^max p^^max _ 

F = - l \ cos^ ri dp dri 

tt J n J n 


Integrating over 3 and expressing 3 in terms of ri. 


BmaxM = cos 1 


rcos il ma x ] 


cos ri 


Since the detector rotates through it radians while 3 goes from 0 to 
^max^)' tiie shape factor can be expressed as 


= if 


’max 


COS^ n cos 


-1 


TT Z 


cos r1 Inax^ 


dn 


cos n 


where the unrestricted view angle 


'max 


= sm 


-1 


+ hQ 
R + h , 


(A2) 


The integral (eq. (A2) ) was evaluated numerically over the range of r) max which 
corresponded to the variation in orbit altitude. The value of D ma x is 55° at 
the mean satellite altitude of 1450 km. 
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APPENDIX 


Direct Solar Radiation 

The solar flux at the Earth can be expressed as 


H 


s 



where H s is the solar flux density at the mean distance of the Sun from the 
Earth and L is the actual distance normalized by the mean distance. 

Now, let a flat-plate detector have an absorptance a and an emittance e 
relative to solar radiation. If the plate is rotated through an incremental 
angle d3 about an axis which is inclined to the Sun’s rays by an angle y, an 
increment of irradiance can be expressed as 


dE s = 



sin y cos 3 d3 


During each revolution, the plate is exposed to sunlight half the time, and the 
irradiance is given by 


or 


1 \ 2 


5 s = H sM ( a / e )s 


sin y p 71 / 2 
\ cos 3 d3 


2tt 


-tt/2 


h \2 sin y 

* ».(-) <“^3 — 


The factor sin y /it can be thought of as a geometric shape factor for a rotat- 
ing flat plate exposed to direct solar radiation. 
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TABLE I.- MONTHLY AND SEASONAL AVERAGES OF RADIATION PARAMETERS FOR 


ESSA 7 DATA PERIOD (SEPT. 1968 THROUGH FEB. 1969) 


Period* 

Solar flux 
density. 

Emitted flux 
density. 

Reflected flux 
density. 

Net flux 
density. 

Albedo 


W/m 2 

W/m 2 

W/m 2 

W/m 2 


Sept. : 






NH 

351 .1 

228.9 

100.9 

21 .3 

0.287 

SH 

317.8 

223.1 

97.4 

- 2.7 

.306 

G 

334.4 

226.0 

99.2 

9.2 

.297 

Oct. : 






NH 

288.2 

231 .8 

83.0 

-26.6 

0.288 

SH 

392.3 

229.3 

128.2 

34.8 

.327 

G 

340.2 

230.5 

105.6 

4.1 

.310 

Nov. : 






NH 

236.1 

225.0 

68.5 

-57.4 

0.290 

SH 

455.1 

231 .6 

148.8 

74.7 

.327 

G 

345.6 

228.3 

108.7 

8.6 

.315 

Sept . -Oct .-Nov. : 






NH 

291 .8 

228.6 

84.0 

-20.8 

0.288 

SH 

388.4 

228.0 

124.8 

35.6 

.321 

G 

340.1 

228.3 

104.4 

7.4 

.307 

Dec. : 



j 



NH 

211 .9 

220.6 

59.5 

-68.2 

0.281 

SH 

486.0 

228.3 

152.7 

105.0 

.314 

G 

349.0 

224.5 

.. 

106.1 

18.4 

.304 

Jan. : 






NH 

224.0 

213.0 

62.6 

-51 .6 

0.279 

SH 

473.7 

228.3 

144.6 

100.8 

.305 

G 

348.9 

220.7 

103.6 

24.6 

.297 

Feb. : 






NH 

268.1 

216.0 

75.0 

-22.9 

0.280 

SH 

425.6 

229.2 

126.8 

69.6 

.298 

G 

346.9 

222.6 

100.9 

23.4 

.291 

Dec .-Jan .-Feb. : 






NH 

234.7 

216.5 

65.7 

-47.5 

0.280 

SH 

461 .8 

228.6 

141 .4 

91 .8 

.306 

G 

348.2 

222.6 

103.5 

22.1 

.297 


*NH, Northern Hemisphere; SH, Southern Hemisphere; and G, Global. 
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TABLE II.- COMPARISON OF ESSA 7 GLOBAL RESULTS WITH OTHER DATA 


Data period 


Average on global scale 


Satellite orbit 


Emitted flux Net flux 

density, density. Albedo 

W/m 2 W/m 2 


Present analysis ESSA 7 

3-p.m. 1 SS 2 orbit 


Sept., Oct., Nov. 


Dec., Jan., Feb. 


7.4 0.307 


22.1 0.297 


Reference 1 


Several satellites , 
different orbits 


Sept. , Oct . , Nov. 


14.0 0.280 


Dec . , Jan . , Feb. 


1 960-61 

1961- 62 

1962- 63 

1963- 64 

1964- 65 


20.9 0.310 


Reference 1 5 


Several satellites, 
different orbits 


Sept., Oct., Nov. 


7.0 0.280 


Reference 8 


Several satellites, 
different orbits 


Dec., Jan., Feb. 


Sept., Oct., Nov. 


1963- 64 

1964- 65 

1965- 66 


7.0 0.310 


6.3 0.301 


Reference 7 


Reference 6 


Nimbus 3, 

12-noon SS orbit 


NOAA 1 , 

3-p.m. SS orbit 


Dec . , Jan . , Feb. 


Oct. 3-17 


Jan. 21 - Feb. 3 


Feb. 18 - Mar. 16 


1964-65 

1968-69 


10.4 0.311 


4.2 0.282 


16.0 0.283 


-12.1 0.351 


'Local time for ascent crossing at 
2 SS , Sun synchronous. 


equator . 
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Figure 4.- Values of shortwave detector response factor (a/e) r for plate 

radiometers on ESSA 7. 
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Figure 6.- Concluded. 
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Figure 9.- Monthly averages of zonal albedo for each month analyzed. 
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Figure 10.- Concluded. 
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Figure IT.- Monthly averages of zonal net flux density for each month analyzed 
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Figure 13.- Global distribution of 

(ESSA 7 data) 
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Figure 14.- Concluded. 
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Figure 15.- Concluded. 
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Figure 16.- Concluded 
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Figure 18.- Geometry for rotating flat-plate detector above Earth. 
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